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ABSTRACT 


Heat  checking  phenomenon  In  mechanical  seals  resulting  from  the 
passages  of  a  single  asperity  over  the  face  of  a  seal  ring  has  been 
solved  with  a  thermomechanical  model  of  a  semi-infinite  medium  sub¬ 
jected  to  the  frictional  heat  source  over  the  moving  contact  zone.  An 
earlier  analysis  of  such  a  problem  was  based  upon  a  two-dimensional 
model  of  a  single  moving  asperity.  Such  a  load  resulted  in  fracture 
which  Is  Initiated  beneath  the  surface  at  the  trailing  edge  of  the 
moving  load.  For  a  better  description  of  the  physical  problem,  a 
three-dimensional  model  Is  examined  In  the  present  study.  A  single 
moving  asperity,  with  a  circular  or  a  rectangular  geometry  of  the  con¬ 
tact  area.  Is  assumed  on  the  mating  part  of  the  seal  face.  The  general 
solution  Is  obtained  for  any  distribution  of  contact  loads.  Numerical 
solutions  are  obtained  for  uniform  and  parabolic  distributions  of  pres¬ 
sure  over  the  circular  and  rectangular  contact  zones.  The  traverse 
speed  of  the  friction  load  Is  considered  the  same  as  the  rubbing  speed 
for  the  current  problem.  The  roles  of  material  properties  and  oper¬ 
ating  variables  are  delineated  In  terms  of  dimensionless  parameters. 

The  analytical  solutions  of  the  mechanical  and  the  thermal  stress 
fields  are  solved  by  means  of  the  double  Fourier  transforms.  The  tem¬ 
perature  field  Is  obtained  by  the  use  of  the  Green's  functions.  Numer¬ 
ical  results  of  the  corresponding  Integral  solutions  are  represented  in 
graphs. 
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Nomenclature 


i  Asperity  characteristic  dimension,  the  half  width  of  the 

rectangular  contact  area  In  the  direction  of  traverse  or 
radius  of  the  circular  contact  area 

>  Half  length  of  the  rectangular  contact  area  perpendicular  to 

the  direction  of  traverse 

:  specific  heat 

)  Differential  operator  with  respect  to  X3 

c  Thermal  conductivity 

Dllatatlonal  speed  ratio  [*V(p/(a  +  2u))^2] 

A2  Shear  speed  ratio  [»V(p/m)^2] 

)0  Average  pressure  over  the  contact  area 

q  Heat  flux  through  the  contact  area 

qQ  Average  heat  flux  through  the  contact  area  (*PfP0V) 

7-j  Traction  over  the  contact  area 

t  Aspect  ratio  (b/a)  or  time 

r  Temperature  field 

{Ul},  u.  Displacement  field 

Y  Traverse  speed  of  asperity  (-xi  direction) 

{ x-f }  Coordinates  fixed  to  the  moving  asperity 

f" )  Fourier  transform  of  a  variable 

Partial  derivative  with  respect  to  x^  coordinate 
»  Coefficient  of  thermal  expansion 

Kronecker  delta 

t  Dimensionless  temperature  field  (»Tk/q0a) 

*  Thermal  dlffuslvlty 


Lam€  coefficient 
u  Lamg  coefficient,  modulus  of  rigidity 

Coulomb  coefficient  of  friction 
9  Mass  density 

Stress  field 

Dimensionless  stress  field 
Mechanical  stress  field 
Thermal  stress  field 

{C,n,c}  Dimensionless  coordinates  (*  x</a) 
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1.0  INTRODUCTION 

The  present  Investigation  addresses  one  specific  failure  mechanism 
In  marine  seals,  that  severely  affects  Its  service  life.  The  seal  Is 
designed  to  locate  along  the  propeller  shaft  to  prevent  leakage  of  sea 
water  through  the  shaft  tunnel.  Figure  1  illustrates  one  assembly  of 
such  device.  In  which  the  seal  consists  essentially  of  two  mating  annu¬ 
lar  rings,  one  of  which  Is  fixed  to  the  tunnel  housing,  another  Is 
mounted  under  spring  pressure  to  rotate  with  the  shaft.  The  mating 
surfaces  of  the  rings  therefore  are  pressed  against  each  other  and 
rubbing  at  a  high  speed.  The  design  nominal  pressure  between  the  seal 
rings  depends  on  the  type  of  vessel  and  the  location  of  the  seal.  For 
instance,  a  low  pressure  of  about  100  kPa  (14.5  psi)  is  sufficient  for 
most  surface  vessels;  for  submarines,  the  pressure  required  could 
Increase  by  orders  of  magnitude.  Had  the  pressure  been  evenly  distri¬ 
buted  according  to  design,  the  life  of  the  seal  will  no  longer  be  a 
serious  problem  even  at  a  high  rubbing  speed.  However,  it  Is  well- 
known  that  the  actual  contact  area  may  only  be  a  small  fraction  of  the 
nominal  area  at  the  design  interface.  In  other  words,  a  low  nominal 
design  pressure  may  very  well  result  In  a  very  high  interfacial  pres¬ 
sure,  thus  a  very  high  dry  frictional  force  In  the  actual  contact 
area.  The  high  friction  would  cause  locally  an  extremely  high  tempera¬ 
ture,  called  "flash  temperature"  by  Archard  [1],  The  local  contact 
area  Is  therefore  otherwise  called  the  "hot  spot"  [2J.  In  severe  cases 
the  temperature  can  be  extremely  high,  leading  to  cracking  of  the  sur¬ 
face  [3].  The  frictional  cracking,  or  heat  checking,  has  been  observed 
to  occur  In  seal  rings  [4].  If  we  use  the  figure  of  10"3  as  area 
ratio  (contact  area/nominal  area)— Burton  [5]  considered  10*  4  as  a 
possible  area  ratio,  a  low  pressure  of  240  kPa  (35  psi)  could  result  in 
a  240  MPa  (35000  psi)  local  pressure  in  the  contact  zone,  a  pressure 
well  within  the  range  of  fracture  Initiation  with  unfavorable  Coulomb 
coefficients  [6], 

The  cause  of  the  localization  of  the  contact  area  will  not  be  con¬ 
sidered  here;  nor  will  the  metallurgical  change  In  the  surface  layer. 


Figure  1.  Mechanical  face  seal. 
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It  will  suffice  to  assume  that  the  localization  of  contact  areas  is  due 
to  some  form  of  asperities,  which  maty  be  fixed  to  any  of  the  mating 
surfaces  or  may  precess  with  respect  to  both.  Kennedy  et  al.  [7], 
using  carbon  ring  against  a  metallic  mating  ring  of  440  Cs.s,  beryllium 
copper  or  52/00  bearing  steel,  showed  existence  of  spot  asperities  on 
the  metallic  rings.  Burton  [8]  reported  on  the  use  of  aluminum  ring  on 
glass  disk,  showing  hot  spot  precesslng  at  a  much  slower  speed  compared 
to  the  rubbing  speed.  Marscher  [9,10]  treated  the  thermomechanical 
stress  state  In  the  part  that  contains  the  asperity.  The  present 
investigation  treats  the  thermomechanical  stress  state  in  the  part  that 
mates  with  the  one  containing  the  asperity  and  the  possibility  of  heat 
checking  therein.  The  formulation  also  allows  extension  to  problems  of 
different  rubbing  speed  and  precesslng  speed  of  the  asperity. 

Because  of  the  relative  size  of  the  seal  material  and  the  moving 
asperity,  the  seal  will  be  represented  by  a  half  space.  The  asperity 
contacts  with  the  material  at  the  otherwise  traction-free  surface  in  a 
small  rectangular  or  circular  region,  which  traverses  the  surface  at  a 
speed  (V)  of  the  order  of  10  m/s  (■’400  ips).  Figure  2.  The  coordin¬ 
ates  {x-j}  are  fixed  to  the  moving  asperity  such  that  xi  points  toward 
the  trailing  direction  of  the  motion,  X2  Is  perpendicular  to  the  tra¬ 
versing  direction,  and  X3  Is  a  depth  measure  pointing  Inward  from  the 
surface  Into  the  material.  Figure  2.  For  the  rectangular  contact  area, 
the  aspect  ratio,  t(«  b/a).  Is  a  parameter  for  three-dimensional 
effect.  For  large  aspect  ratio,  the  load  Is  effectively  a  moving  line 
load  that  has  been  represented  by  a  two-dimensional  plane  strain  solu¬ 
tion  [6].  In  the  present  report  the  pressure  and  tangent  force  are 
arbitrarily  distributed  over  the  contact  area  for  the  general  solutions 
In  the  Integral  form.  Numerical  solutions  will  be  obtained  specifical¬ 
ly  for  the  uniform  and  the  symmetrical  parabolic  distributions. 

Since  the  speed  of  traverse  Is  much  smaller  than  the  Rayleigh  wave 
speed,  which  for  a  steel  based  seal  material  is  approximately  2800  m/s 
(~11  x  10**  Ips),  no  wave  propagation  phenomenon  Is  considered.  How¬ 
ever,  for  the  same  material  and  a  half  contact  width  of  the  order  of 


(a)  Rectangular 
Contact  Area 


(b)  Circular 
Contact  Area 


1  mm,  the  Pecl£t  number  (Pe  *  Va/ie)  Is  of  the  order  of  500,  of  which 
the  magnitude  Is  quite  sufficient  to  necessitate  the  consideration  of 
the  dynamic  effect  of  the  traverse  speed.  It  is  assumed  that  the 
material  remains  predominantly  elastic  and  hence  that  the  uncoupled 
theory  of  thermoelasticity  holds.  Later  it  will  be  verified  that  the 
maximum  principal  stress  will  be  under  the  yield  stress  except  at  a 
very  small  neighborhood  where  the  ultimate  strength  is  reached.  The 
postulations  allow  us  to  treat  the  mechanical  stress  state  and  the 
thermal  stress  state  separately.  The  combined  effect  will  then  deter¬ 
mine  the  possibility  of  fracture  initiation.  The  general  solutions  for 
the  mechanical  stress  state,  the  temperature  field  and  the  thermal 
stress  state  are  derived  with  the  Fourier  transform  method  and  the 
Green's  function  method.  The  numerical  data  are  computed  with  a  com¬ 
bination  of  analytical  integration  and  numerical  method  to  avoid  heavy 
computation  time  and  stability  problems. 


2.0  STRESSES  FROM  MECHANICAL  EXCITATION  {a7  } 


_ L* _ 

_Bds1c_^overni_n3_i3Mil5i2IL5.*  The  working  equations  come  from  Cauchy's 
law  and  Hooke's  law.  In  terms  of  the  moving  convective  coordinates 
{xj}  and  with  the  absence  of  body  forces,  the  acceleration  in 
Cauchy's  law  will  have  only  the  convective  terms.  Hence,  we  have 


3. 
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Vi 


°ij  ■  *Wlj  +  1,(3i“j  *  Vl*  ’  l.j.k  -  1.2.3  (2 

where  u^  is  the  displacement  field,  o1^  is  the  mechanical  stress  field 
3^  denotes  partial  derivative  with  respect  to  ,  6^.  is  the  Kronecker 
delta,  V  is  the  asperity  traverse  speed,  (p,  x,  u)  are  material  con¬ 
stants.  The  method  to  solve  (1)  and  (2)  will  be  similar  to  the  one 
used  by  Eason  [11]. 

Equations  (1)  and  (2)  are  conveniently  solved  by  the  method  of 
double  Fourier  transform  with  respect  to  (xi,  X2) 


T(7lfx2,x3) 


ifox^+x-x,,) 

f(x1,x2,x3)  e 


dx^  dx2 


(3 


With  the  use  of  the  relation 


a"  f  *  (-1  *r)n  T  ,  r  •  1,2  , 
we  get 

-17r  °jr  +  0  °J3  *  -pv2  *1  “j  >  J  -  1-2-3  ;  r  *  1>2  >  <4 


1Summation  convention  is  used  for  repeated  indices  of  roman  minuscules 


o  a  -IX  x.  u,  6  -  1  u(x  u  +  x  u  )  +  XDu^fi 

rs  k  k  rs  v  r  s  As“r'  3  r 


°3r  "  *,uxru3  *  uDur 


k,r,s  -  1,2  ,  ( 


-Ux-u.,  +  (X  +  2w)  Du, 


subject  to  the  boundary  conditions 


°31  *  -*1  at  x3  ’  0 


Uj  ♦  0  is  x,  ♦  •  ,  ( 

where  the  transformed  quantities  are  denoted  by  a  superposed  (“ ). 
is  the  traction  on  the  surface  X3  *  0,  D  denotes  derivative  with 
respect  to  X3  and  the  Index  1  *  1,2,3. 

Displacement  and  Stress  Fields.  Solution  of  the  set  of  ODE  (4,5)  fol 
lows  by  combining  (4)  and  (5)  first,  resulting  in 

[D2  +  (M2  -  e2)  T2  -  x2]  ux  -  xjx2(s2  -  l)  u2  -  ixx(e2  -  l)  Du3  =  o 

*2 ( ®2  *  1)  “i  +  CD^  +  (M^  -  1)  x^  -  32x|]  U2  -  172(32  -  1)  DlTj  =  0 

-ix^P2  -  1)  Ouj  -  ix2(02  -  1)  Du2  +  [32D2  +  (M2  -  1)  x^  -  x2]  u3  a  0 


in  which 


M1  *  V(rT7il)  *  m2  *  v(^}  ;  •  *  (*  +  2w)/u  . 


Equation  (8)  Is  a  system  of  three  simultaneous  homogeneous  linear 
equations  for  three  unknowns  {IT| } .  For  nontrivial  solutions  the 
determinant  of  coefficients  should  be  zero  Identically,  that  Is 


D2  +  (M2-B2)  x2  -  x\ 


-  x1x2(B2-1) 


-x^B2-!) 


■Ix^B2-!)  D 


1x2(B2-1)  D 


-Ix^B2-!)  D 


D2  +  (m|-1)  Xj  -  B2  x| 


■1x2(B2-1)  0 


B2D2  +  (m|-i) 


*  B2(D2  -  n2)  (02  -  n2)2  *  0  , 


where 


nr  *  +  »  r  *  1,2  . 

Hence,  Equation  (8)  results  In 

(0  -  nj)  (0  -  n2)  (u j ,  u2>  u3)  ~  0 

whose  solutions,  with  the  consideration  of  the  regularity  condition 
(7),  are 

-n  i  x-j  -n9x-> 

U1  ’  Ale  i  +  (B1  +  cix3^  e 

-n,Xo  -n9x, 

u2  »  A2e  +  (B2  +  C2x3)  e 


u3  *  A3e  °lX3  +  (B3  +  C3x3)  e'n2x3 


t 


( 


where  Ai,  A2,  As,  Bi,  82,  83,  Ci.  C2,  C3  are  independent  of  X3. 
Substituting  Equations  (10)  into  (8)  gives  certain  relations  between 
these  quantities  as  follows. 


A 


r 


r  -  1.2  , 


X«B«  +  Xnfi; 

b3  ■  • 


and 


Cj  *  0  ,  j  *  1.2,3  . 


Thus.  Equation  (10)  can  be  rewritten  as 


__  1x^3  ”^1*3  ”^2*3 

U1  "  ~H7 —  e  +  Bj  e  , 


<*2*3  -"i*3  .  „  _-n2x3 

y  3  c  * 


2  n, 


B2e 


-n,  Xo  t  HlpXo 

U,  ■  A,e  1  J  -  —  (x^Bj  +  x2B2)  e 


'3  3 


(ID 


Expressions  for  the  transformed  stress  components  may  be  obtained 
by  the  substitution  of  Equations  (11)  into  (5).  In  particular,  we  have 


XMJxJ  -  2vnJ  ^ 
nj  *3 


-n1  x<%  "IIaX^ 

e  1  J  +  ^(xjBj  +  72B2)  e  £  J 


o32  »  -^ux^e  1  3  -  UjXgBj  +  (n|  +  x|)  B2]  e  23 


°31  3  -^MXjAj  e  1  3  -  —•  [(n|  +  Xj)  B1  +  x^Bg]  e  2  3  .  (12 

In  view  of  the  boundary  condition  (7),  A3  and  Bi,  B2  can  be 
readily  solved  and  expressed  in  terms  of  7j. 

2mFA3  *n1[(x^  +  x|--j  M^)  73  +  IngCx^  +  7272)]  , 


2uFBj  »  -fx^j^n^  +  i-  {[xgfx^  +  “  7  m|xj  -  2njn2)  +  n2(x|  +  7? 


”  7  ^*1^  ”  xlx2^*l  +  *2  ”  7  ^2X1  ”  ^2^  * 


2uFB2  »  -ix^n^  +  -p--  {-XjX2(x^  +  x|  -  M2Xj  -  2n^n2)  TT^ 


+  [x^(x^  +  x|  -  MgX^  -  2n^n2)  +  n2(x^  +  7^  -  M27^)]  1^}  . 


Following  the  substitution  of  Equations  (13)  into  (11),  we  obtain 

_  o  _  -n,x3  (  _ 

2uF  *  (lx ^2^3  *  Xxn2^1  *  x 2^ 2^2^  ^  *  |1x 2^3 

r  *2  "1  *i  X9  )  _n?X3 

-  — —  (n,  -  2n-.ru)  +  n«n«  7,  +  -  (n-j  -  2n,n0)  7-  /  e 


2yF  u *  (IXgnjUj  •  x^x^n^R*^  •  x^n^R^)  ® 


xix? 

+  -B^<"3 


f*l 

-  2n1"2)  *1  -  Tq 


jlx^njnj 


(n3  -  2^2)  +  n2n3  /  e 


2m F  u3  *  (nin3^3  +  Ix^n^!  +  IXgnjngR^)  e 


where 


-  Cnl Cxf  +  *f )  ^3  +  ^xln3^l  +  ^x2n3^2^  e  2  3  » 

1 

n3**l  +  *2"7^2*l  *  F  =  n3  -  nin2(xl  +  *f)  * 


Stress  components  can  then  be  determined  by  the  substitution  of  (14) 
Into  (5)  In  accordance  with  the  new  variables  as  follows 

N1  “  n3  ”  2nin2  »  N2  *  Ni/n2  »  ^  *  n3^3  +  in2^xl^l  +  "*$2}  * 


thus  we  have 


2uF  *  TT  [(x^  +  x|  *  nj)  *  +  2ux^]  e  1  3  -  2M{x^n1n2Tr3 
+  Ixj[x|n2  +  ngn3]  1^  -  Ix^xyigR’g}  e  2  3  , 


_n2x3 


2mF  0^2  a  *  C(xf  +  xf  -  h2)  *  +  2mx|]  e’"1*3  -  2p{^n1n1lT; 


*  lx ^  1  x2Lx ?N 2  ^  ^2^3^  ^2^  ® 


’n2x3 


1  0 


”n2x3 


2uF  0^2  *  b[(x^  +  x£  -  n£)  A  -  2un£j  e 

+  2un2[(7^  +  72)  +  in^fx^  +  ^"2^2 ) ^  e 

2wF  oJ2  *  2u7172'ff  e  1  3  -  2y{x^x2n^n273  +  ^  x2Cn2 


~  (xj  -  x2)  N2]  7^  +  ^  x^[n2n3  ^1  "*  x2^  ^2^  "^2^  ® 

2yF  "0^3  =  -2wl72n^¥  e  *  3  +  u{l72n^(x^  +  x|  +  n2)  7^ 

-  2x^7^ -  [7^  +  n3(n2  +  *f)]  e  2  3  * 


2uF  0^  *  -2111x^7  e  1  3  +  u{ixjn^(7j  +  x|  +  n2)  73 

-  [x|Nj  +  n3(n2  +  7^)]  7^  -  e  2  3  . 

The  displacement  and  stress  fields  are  to  be  left  In  the  Fourier 
transform  expressions  (14)  and  (15).  The  actual  fields  are  to  be 
obtained  through  Inverse  transform  after  the  substitution  of  the 
boundary  values  {7<}. 


3.0  STRESSES  FROM  THERMAL  INPUT  {o[ .} 


Temperature  Field.  The  heat  equation  with  constant  thermal  properties, 
assuming  quasi -steady  state  and  no  heat  generation  In  the  medium,  as 
expressed  In  terms  of  the  convective  coordinates  {xj},  is 


»iiT-*  V 


where  *  is  the  thermal  diffusivlty. 

The  boundary  conditions  at  x 3  *  0,  are 


!-q(xp  X2)  ,  In  the  contact  region 

(17) 

0  ,  elsewhere 

T,  3^T  ♦  0  as  (x^)1'2  -  -  .  (18) 

Note  that  T  Is  the  temperature  above  the  stress-free  ambient,  k  is 
the  thermal  conductivity,  q  Is  the  heat  flux  (-RiV  *  MfP V),  trfiere 
Uf  Is  Coulomb  coefficient  and  p  Is  normal  pressure  distribution  on 
the  surface  X3  3  0.  It  is  assumed  that  the  mating  seal,  that  contain¬ 
ing  the  asperity.  Is  an  Insulator. 

The  double  Fourier  transform  of  Equations  (16)  through  (18) 
become,  respectively, 

[D2  -  (s2  +  In)]  T  «  0  ,  (19) 


kOT  -  TJ 


-q(xj,  x2)  ,  In  the  contact  region 


elsewhere 


T  ♦  0  as  x,  ♦  •  . 


where  s2  -  x|  +  x|  and  n  =  -(V/ie)  xj.  Now  let  T  *  +  iT2  and  Q/k  * 

Tj  ♦  i^2*  then  Equations  (19)  and  (20)  result  In 

(D2  -  s2)  T1  -  -n?2  , 

(D2  -  s2)  T2  «  nT1  (22) 

and 

°T1  *  ?1  • 

DT2  «  Pg  at  x3  -  0  .  (23) 

Also,  (21)  becomes 

Tp  T2  *  0  as  x3  ♦  -  .  (24) 

Equation  (22)  can  be  expressed  as  a  fourth-order  homogeneous 
ordinary  differential  equation  which  Is 

[D4  -  2s2D2  +  (n2  +  s4)]  Tr  »  0  ,  r  »  1,2  .  (25) 

In  view  of  the  regularity  condition  (24),  the  solution  of  Equation  (25) 
has  the  form 

Tr  ■  e"“z[Af  cos  0z  +  Bp  sin  0z]  (26) 

where 


14 


Following  the  substitution  of  the  expression  (30)  into  (26)  and 
according  to  T  ■  Ti  +  iT2»  we  obtain  T  as 


_  -»x« 


T  *  e  [(C,  cos  8x3  +  C2  sin  0x,)  P 


3'  rl 


+  (C3  cos  0x3  +  C4  sin  8x3)  P2]  , 


where 


-2«  +  in 


2u(uf  +  8^) 


,  29‘  +  In 

*9  5 - T~  » 

c  29 (or  +  8^) 


-n  -  12^ 
2»(m2  +  8^) 


-n  +  i 28* 


4  29  (tf  +  0  j 

Again,  the  temperature  field  Is  left  in  the  Fourier  transformed  expres 
sion,  which  will  be  used  in  solving  the  thermal  stress  field.  To  be 
complete  and  also  for  comparison,  the  Green's  function  solution 
approach  for  the  temperature  field  is  given  in  the  following. 

The  Green's  function  solution  can  be  expressed  as 


g(xj  -  xj,  x2  -  x£,  x3,  t  -  t') 


4pC[*ic(t  -  t')]; 


C(X,  -  xi)  -  V(t  -  t')]2  +  (x,  -  xA)2  +  x2 


which  is  the  temperature  at  time  t  at  the  location  (xi,  X2,  X3)  due  to 
the  unit  heat  flux  Input  emitted  at  time  t'  at  the  location  (xj,  x£, 
0).  It  can  be  easily  shown  that  Equation  (32)  satisfies  the  governing 
heat  equation  (16). 

The  temperature  field  can  be  obtained  by  the  convolution  integral 


^(xl*x2*x3»t)  =  /  q(xlsx^)  Kq(x^  -  x^, 
J  B>> 


Xg  -  x£,  x^,  t  -  t')  dxj^  dx£ 


where  Indicates  the  heat-flux  Input  region  and 


'll  ^  '  xi' 


x2  *  Xj,  t  -  t')  dt'  . 


Substituting  (32)  Into  (34)  and  changing  the  variable  by  letting 
t  ■  R/2/ie(t  -  t').  Kg  becomes 


V(Xj-Xj)/2x  rm  2  1  9  70 

K  •  9  „,-■■■■  I  e-t2-(V2R2/16«2,2)  dt 

0  *3/z  kR  JR/2/n 


where  p 


density 
specific  heat 


C(xx  -  xj)2  +  (x2  -  x£)2  +  x2]1/2 


Considering  the  steady-state  solution,  i.e.,  t  ♦  «,  Equation  (35) 
results  In 


Accordingly,  the  steady-state  temperature  field  can  be  obtained  as 


KXi.X2.X3)  q(xi»x^) 

yB  H 


-V[R-(xrxp]/2»c 


dxj  dx£  (37) 


which  will  be  employed  to  give  the  numerical  results  for  different 
cases  of  application. 


Thermal  Stress  Field.  The  governing  equations  for  the  thermal  stress 
solutions  are  based  on  the  theory  of  quasi-static  uncoupled  thermoelas¬ 
ticity.  The  governing  equations  are  then 


t»V2Uj  +  (*  +  v)  •  u)  *  (3X  +  2u)  a3jT 


(38) 


and 


'ij  *  XWlj  +  w<3ju1  +  31uj} 


(3*  +  2u)  aT6 

*  v 


(39) 


subject  to  the  boundary  conditions 


'31 


0  ,  1  *  1,2,3  at  x^  *  0  , 


(40) 


ojj  ,  ui  ♦  0  as  (xixi)1/2  ♦  -  . 


(41) 


The  double  Fourier  transforms  of  Equations  (38)  through  (41)  are 


[D  -  UjX^  +  Xg) J  -  ia2x^Du3 


_  _ _  p  _ p  _ p  _  _  _ 

-  a2XjX2Uj[D  -  (ap^  +  xp]  u2  -  ia2x20u3  3 


_  _  __  22—  — 
-  ia^jDu^  -  ia2x20u2  +  (a^D  -  s  )  u3  3  a^OT 

o},  __  __  _  - 

-ii  3  -i(a1x1u1  +  rx2u2)  +  rDu3  -  a3T  , 


*L  __  __  _  - 

*  -'(rxjUj  +  ajx2u2)  +  r0u3  -  SjT 


3  -ir(xjU^  +  x2u2)  +  a^Du3  -  a3T  , 


-r  ■  -i(x2ui  *  xiu2>  > 


"T  *  Du2  -  1x2u3  • 


^  _  __ 

—  ■  0ul  -  1xlu3  • 


°31  ’  0  • 


i  3  1,2,3  at  x,  3  0  , 


?j.  ,  u1  ♦  0  as  (Vi)1/2  ♦  -  , 


where  r  3  X/y,  a  3  r  +  2,  a„  3  r  +  1  and  a,  3  (3r 


Equations  (42)  through  (44)  form  a  system  of  three  simultaneous 
linear  differential  equations.  The  dependent  variables  11,  are  to  be 
solved.  The  system  can  be  rewritten  as 

Au.  *  4^  ,  i  *  1*2,3  ( 

where  A  is  a  differential  operator  that  is  the  determinant  of  the 
system. 


D2  -  (a^j  +  x|) 

-Vlx2 

-ia2X jD 

a  = 

- *2*1*2 

D2  -  (ajxf  ♦  xf) 

-iagXgD 

-ia^D 

-ia^D 

3l02  -  s2 

-  ajID2  -  s2)3  ( 

and  Aj  's  are 


”i a3x 1T 

•Wz 

-ia2X^D 

41* 

"*a3x2^ 

D2  -  (SjX2  ♦  X2) 

-ia^D 

a30? 

-ia^D 

axD2  -  : 

1  D2  -  (a,xf 

+  x?)  -ia^.T 

-ia,x\l 

(a^  +  x|) 

'a2*lx2 

"ia3xiT 

-a 

02  -  U,x§  *  x2) 

"^a3x2^ 

-ia^D 

-ia^^ 

a30T 

The  solutions  obtained  through  Equation  (53)  are  particular  solutions 
which  are 

_ p  a-X.  -«x, 

>  *  (Gr  cos  0X3  +  Hr  sin  0X3^  e  »  r  s  1*2  * 

d  a,  -ux, 

“3  -  -Zip*  <G3  cos  9x3  +  H3  s1n  9x3>  e  (S 

where 

G1  3  G2  *  ®2  “  *  Hj  *  Hg  3  -Ag  +  1A^  , 

63  =  -UuGj  -  QHj)  ,  H3  *  -IfwHj  -  ©Gj)  . 

Ai,  A2,  Bi  and  82  are  shown  in  Equation  (30).  The  complementary 
solutions  can  be  determined  by  =  0.  In  view  of  Equation  (54)  and 
considering  the  regularity  condition,  we  obtain 

-  X,  -SX. 

«  -i  (hj  +  h^)  e  , 
x2 


-sx3 


u?  3  (h1  +  h9x,)  e 


t 


u^  ■  — —  C a2 s ( h ^  +  h2x3)  +  (a2  +  2)  h2J  e  .  (56) 

*2*2 

With  regard  to  Equations  (55)  and  (56),  the  general  solutions  are 

u^  *  u^  +  "u^  ,  i  =  1,2,3  .  (57) 

Two  unknowns  hi  and  h2  will  be  determined  by  substituting  Equation  (57) 
into  Equations  (47)  and  (49)  and  employing  the  boundary  condition  (51), 
thus  we  have 

•j  ^ 

hj  -  — ^2 —  [a^se  Hj  +  (s2  -  a^w  +  a^2)  Gj]  , 

Ea^agS  <*>0 

h2  *  2a^6  C"s8  H1  "  "  s“  +  aie2)  Gil  •  (58) 


Following  the  substitution  of  the  displacement  solution,  (56) 
through  (58),  Into  Equations  (45)  through  (50),  the  thermal  stress 
field  can  be  readily  obtained. 


•0>x7  -sx, 

iK-bjHj  +  b^)  e  3  *  (-b3H1  +  b4Gx)  e 


-T 


°99  r  _UX?  "SX3 

-f- »  1  (-bgHj  +  bgGj)  e  3  +  (-b?H1  +  b^)  e  3  , 


oT 


P  >  if(-b0H,  +  blnG,)  e"”*3  +  (b-.H.  -  b„G.)  e  S*31 


* 


—  *  i^-b13Hl  +  b14Sr  e  +  ^’b15Hl  +  b16Gl ^  xlx2  e 

L  r  -wx,  -sx^ 

F*  “x2  [<b17Hl  -  b18Gl>  e  3  -  (bi9Hi  -  b20Gi)  e 

3i  _  r  -««*  “sxi 

u  a  _Xl[(b17Hl  -  b18Gl>  e  “  (b19Hl  '  b20Sl>  e 


-(n  +  i2u2}  F  +  (2<*>2  -  In)  F 

i  * - =5 - 5 - -  *  i(C,F,  +  C/,)  , 

1  2u>(oid  +  QZ)  1  1  3  Z 


(n  -  1202)  +  (292  +  in) 

20(w2  +  9Z) 


*  “  1  ( C 2^ 2  +  C4F. ^)  * 


1  "  ajlff  <2“®  cos  8x3  -  *1  s,n  9x3> 


2  ■  ajis  <*!  cos  9x3  *  2“°  s1n  0x3>  • 

3  u  8282  a2  /  * 


*16  "  aJOT  [  a2s  +  (s  "  u)  x3j  • 
’17  "  ap-  <-6  cos  9x3  *  “  s,n  9x3> 
>18  '  ap  <“  cos  9x3  +  9  s1"  9x3>  • 


919  *  ipr  t'1  +  sx3>  • 


ip  t-“  +  s(“  -  s>  x33  • 


Note  that  7*  and  T2  are  the  real  and  imaginary  part  of  IJ/k 
Therefore 

’’l  ■  2il T  /„  «xl»  x2>  “s  (Vl  +  X2X2)  dxl  dx2  • 

jbh 


TT  _  * 

v2  l 


/  Q(x1, 

JQ..  l 


Xg)  sin  (xfjXj  +  X2X2)  dx^  dx2 


For  the  case  Tz  »  0,  Equations  (59)  become 


fL  -«x7  -sx, 

"y~  *  7^  (b^C^  +  ^2^2^  ®  +  ^3^1  +  ^4^2^  ® 


°L  r  -ux,  -sx ,  1 

-rr-  -  7,  (t>eC,  ♦  b<C,}  e  3  *  (b,C,  ♦  b„C„)  e  J  , 


4.0  NUMERICAL  RESULTS 

Previous  sections  have  shown  the  general  expressions  for  the 
transformed  mechanical  and  thermal  stress  fields.  Applying  the  Inverse 
Fourier  transform  and  following  through  the  simplification  of  the 
multiple  integrals,  results  for  several  cases  can  be  obtained 
numerically. 

4.1  A  Moving  Rectangular  Contact  Area  of  Uniform  Pressure  (Case  1) 
Mechanical  Stress  Field 

The  tractions  of  the  mechanical  loading  are 

°32  '  -R2  *  0  •  °31  *  -R1<V  x2*  °> 


xj  <a,  x2  <b 
elsewhere 

R1  *  -“^3  • 

The  corresponding  transformed  boundary  conditions  can  be  readily 
obtai ned 


°33  3  "R3^xi»  x2*  ' 

where 

!P0  «  constant  , 

0  , 


x2* 


1  .,(7lxltx2x2>  ,  d, 

tw  J  m  R3  e  dxl  dx2 

y /;  /;  *.  *, 

i  [  /.: -,j  [  /:  •%i  *, 


(Equation  continued  on  next  page.) 


(6 


sln(x^a)  sin(x2b)  , 


*xlx2 


*  TTg  *  0  . 


Following  the  substitution  of  (62)  Into  (15)  and  applying  the  inverse 
transform,  the  mechanical  stress  components  can  be  expressed  as 


■  Kji  -  Kj) 


where  the  first  term  on  the  right-hand  side  of  this  equation  denotes 
the  stress  components  resulting  from  the  normal  loading  ^3,  and  the 
second  the  frictional  loading  ffi.  The  expressions  are 


n  po  f2'  cot  ♦  /""  f 

’ll  7 Jo  -^-Jo  [ 


r3(l  +  \  n| 


9 

Mf)  e  1  6  -  y1y2  e 


-sy2x3 


•  W(s,$)  ds  d$ 


n  s  Po 

22  7 


J*  TT  'sin  co's  [t3<i  82  -  1)  Mf  cos2  ♦  ♦  sir,2  ♦] 


-sy,x3 
e  -  yxy2  sinfc  4>  e 


2  "sy2x3 


W(s,t)  ds  d* 


_n  _  o 
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i t  h 


YiY2  e 


'SY2x3 


W(s,*)  ds  d* 


,n  .  p0  [*  i  /•-  r 

*12  "7  L  TT  L  r 


-SYix3 

3  e  *  y1t2  e 


sy2x3 


W(s,*)  ds  d$  , 


cos  ♦jo  r 


j  W(s,4>)  ds 


//[• 


"sy1x3  "SY2x3l 
e  J  -  e  £  J  W(s,*)  ds  d+  , 


7 Jo  |t2  cos2  ♦(!  +  -  «^)  e' SV3 

“  2YlV  2  -SY?x3 

- sin  ♦  +  y2y3  e  >  W(s,4>)  ds  d<(> 

t 

f  |y2^7  “  1)  Mj  cos2  ♦  +  sin2  ♦]  e  Yl 

f(Y3-2Y1Y2)  2  1  -sy-xJ 

+  {--y2  '  sin  ♦  I  e  23  H(s,  +  )  ds  d* 

*fo  |C(*  p2  *  1]  "i  cos2  ♦  •  ^ 


’sy1x3 

e  +  Y3  e 


W(s. ♦)  ds  d$ 


fm  (  2  .  _sy1x3  1  r(Y3  “  2y1y2> 

Vo  r 

(sin2  ♦  -  cos2  ♦)  +  y,1  e  w(s,4>)  ds 


*»  *  -rjo  o  ie 


J  W(s,$)  ds  d$ 


f  Vo  f2"  i  f~  r 
J31  *  -7-  Jo  xmrTc*  1 J0  [1 


ni  'SV3 
V2 cos  ♦ e 


.SY  ^ 

-  (y|  -  y^2  sin2  <j>)  e  23  W(s,4>)  ds  d<> 


where 


W(S,*)  *  sin(sa  cos  »)  sin(sb  sin  <t>)  e'1SY  cos^0  " 

Y  3  (x2  +  x\)l,Z  ,  y3  *  (1  -  J  M2  cos2  ♦)  , 

Yt  *  (1  -  ^  cos2  *)1/2  ,  t  «  1,2  , 


2  -1  x2 
H  =  Y3  -  YjYg  ,  0  =  tan  A  ~  . 


Equations  (63)  can  be  reduced  to  single  integrals  by  using  the 
techniques  of  change  of  variables  and  the  following  equalities: 


s1o(s.acos  sin  jj  e*sV3  cos(sT  cos(9  .  t)) 

1  1  "SYtx3  r 

*  e  {[cos  s(a  cos  4>  -  b  sin  ♦  -  y  cos(e  -  $)) 

+  cos  s(a  cos  $  -  b  sin  +  y  cos(0  -  ♦))] 

-  [cos  s(a  cos  ♦  +  bsin^-Y  cos(0  -  ♦)) 

+  cos  s(a  cos  $  +  b  sin  $  +  y  cos(8  -  ♦))]}  , 


miLSBl  »?sSln(sb  sin  ♦)  #  ’t*3  sin(s  y  cos(9  .  #)) 

1  1  "^Yt*3 

*  e  {[sin  s(a  cos  ♦  -  b  sin  ♦  +  y  cos(e  -  $)) 

-  sin  s(a  cos  ♦  -  b  sin  ♦  -  y  cos(e  -  ♦))] 

-  [sin  s(a  cos  ♦  +  b  sin  ♦  +  y  cos(9  -  ♦)) 

-  sin  s(a  cos  ♦  +  b  sin  ♦  -  y  cos(0  -  $))]}  , 


/  j  e"sd  sin(se)  ds  *  tan-1  (|)  , 

fo  T  e'Sd  C0S(Se)  dS  mfo  s  e'Sd  dS  -  7  ,n  • 

The  term  -i  e“sd  ds  cancels  out  during  simplification.  Thus,  the 
stress  fields  are  expressed  as 

°U  =  -£zf0  Cy3(1  +  7  ’i  ’  Ml)  S1  "  Y1Y2S2]  d*  » 

°22  *  ~Z?  f  TTsTiTZ?  *[y3(7  02  '  X)  "l  cos2  *  +  si°2  *] 


•  Sj  -  yxy2  sin2  <p  S2}  d* 


,n  =  -po  r  i 
33  T?  Jo  H  s1 


CT^S-,  •  y 


H  sin  2?  li331  ‘  t1t2°2 


*12  "  77  L  ¥  ^y3S1  "  Y1Y2S2]  d*  » 


c 


HsTnT  L"1  2J 


(  ' 

'  F1TT7  !t2  cos2  ♦»  +  7  4  -  &  W1 

.  [(T3  -j\V  s1n2  *  +  Vj]  „2J  dt 

*  rsTST  jv<?  ®2  • 11  A  cos2  ♦  *  s1n2  *]  H1 

r (y3  -  2y,Y2)  2 

+  --T  --  --  Sln^  ♦  W2[  (!♦ 

> 


*  Y 


rsIFT  tC(7  82  -  D  ^  “s2  ♦  -  Tl]  “l  +  1,3W2)  d*  • 


*  Y«  (  o  1  r^3  ~  2Yi Y2 ) 

,  fnJsT  “s  ♦  M1  -  7  /  - 


•  (sin2  ♦  -  cos2  ♦)  +  yJ  W2>  d«|» 


*  TiYp 

.  “H  CS1  '  S2]  d*  • 
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■ufPo 

1?" 


f 

Jo 


1 


H  sin"  « 


1 —  [t^2  cos  *  si 


<Y|- 


y1y2 


sin^  <>)  S9]  d<|> 


(64) 

in  which 

{[a  cos  ♦  +  b  sin  <fr  +  y  cos(e  -  4>)]2  +  } 

S,  *  In - 5 - 44- 

1  {[a  cos  ♦  -  b  sin  ♦  -  y  cos(0  -  ♦)]  +  y^x^} 

{[a  cos  $  +  b  sin  $  -  y  cos(e  -  <j>)]2  +  y?x?} 

+  In - n - 44- 

{ [a  cos  4>  -  b  sin  <fr  +  y  cos(0  -  $)]  +  y^x^} 


W  =  tan-1  a  cos  ♦  ~  b  S1n  ♦  *  T  cos(9  -  <t>) 

1  Ytx3 

+  tan-l  a  cos  $  +  b  sin  <)>  -  y  cos(6  -  $) 

Ytx3 

-  tan“^  a  cos  ♦  -  b  sin  <t>  -  y  cos(6  -  ») 

Ytx3 

tan-l  a  cos  <fr  +  b  sin  <fr  +  y  cos(e  -  $) 


t  -  1,2  . 

The  integrals  appearing  in  Equations  (64)  are  suitable  for  numeri¬ 
cal  evaluation  at  any  point  inside  the  solid.  However,  most  integrals 
are  not  defined  at  <fr  =  0,  w/2  or  *.  Limiting  processes  are  imperative 
so  that  singularities  can  be  removed. 
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Thermal  Stress  Field 

Replacing  Q(xi,X2)  in  (60)  by  expressions  (20)  and  (17),  we  obtain 


"P"l  3  /  R(xj»  X2)  cos(xlXl  +  X2X2 )  dx^  dX2 

J  bh 

-ufPV  fa  fb 

*  J  I  cos(x1x  +  x2x)  dx2  dxx 


-MfPQV  sin(x2b)  fa 


i(x2b)  fa 

— = —  /  cos(xjx)  dx^ 

X,  J- a 


-w.  1 


•  _ __  sin(x^a)  sin(x2b) 


72  “  fx  q(xl*  X2J  sin(xlxl  +  X2X2)  dxl  dx2 


*fPnV  /*a  /*b 

^/a  /b  Si"°‘1 


x  +  XgX)  dx2  dxj 


=  0  . 


Substituting  the  expression  7i  into  (61)  and  applying  the  inverse 
Fourier  transform  yields 


0T1  -4UeP  V  /*“  f°°  -<*w, 

-T  *  -^J2-  J0  J0  fl  C0S<X2X2>  [e  (blDl  *  b2D2> 

-SX-,  ]  _ 

+  e  (b,0,  +  b.D,)  dx,  dx. 


Jo  Jo  fl  C0S(72x2)  [e  <b5Dl  +  b6°2] 


+  e  vb7°i  +  b8°2^  ^*1  dxj 


>1,  -4w-PV  f-  /*-  -«x, 

—  *  -  JQ  yo  fi  cos (x^2)  e  (bgD1  +  b1QD2) 


"  e  (bllDl  +  dx^  dx. 


>L  -4u*PnV  /*-  /*-  f  -«x, 

-T  ■  1  J0  fl  Sl"<V2>  [e  '-b13°3  *  b14C 


“  xlx2  e  (bi5°3  +  bl6°4^  ***1  ^*2 


»L  *4w-P  V  /*-  /■-  f  -UX-, 

*■ ■ ^rio  Jo 7*  le  (b 


17°5  +  b18D6} 


®L  -4yfP  V  /*«  /*«  _ 

- a  -?r  Jo  Jo 7ifi  C0S(7*’ 


"e  (b19°5  +  b2(W  ^*1  ^ 


e  (bi7°7  +  bi8°8^ 


-  e  (bi9°7  +  b20D8^  ^*1  ^ 


In  which 


9 


fi  =  -  =-=-  sin(x1a)  sin(x2b) 
xlx2 

Di  ■  *  XT?  [“  C0Spri,Ii)  +  h  s1n(*ixi>]  • 

°2  3  ~T~?  cos(*ixi)  -  75  sin^Xj)]  , 

°3  s  -  ~T~2  sin(*ixi)  +  73  cosOTjXj)]  , 

°4  a  7-*  q2  C9  slnfxjXj)  -  cosfx^)]  , 

05  *  g  1  g  [-  -grr  SlnfXjXj)  +  u  cos  (Six,)]  , 

Dg  s  — 2  1  02  sinCx^x^)  +  0  cos(xjx^)]  , 

°7  *  M7T'e'?  cos(*lxl)  ’  “  s1n(7jXj)]  , 

08  *  J  l  e2  ^  cos(xixi)  +  9  sin (jT^x x ) 3  . 

Due  to  the  complexity  of  the  expressions,  no  existing  automatic  Inte 
gratlon  scheme  for  the  above  Integrals  are  efficient.  Hence  the 
Gausslan-Laguerre  quadrature  formulas  are  used  for  the  numerical 
approximation. 


4.2  A  Movlnq  Disk  of  Uniform  Pressure  (Case  2 


Mechanical  Stress  Field 

The  pressure  distributions  are  assumed  to  be  of  the  form 

x2»  *  °32  *  "R2  *  ®  =  *2*  0) 


where 


?  o  p  o 

PQ  ,  y  *  xj  ♦  Xg  <.  a 


0  ,  elsewhere 


R1  *  -»f*3  * 

The  transformed  boundary  conditions  can  be  obtained  by  using  Equation 
(3). 


ffjfZj,  x2.  o)  - 


r.2.,2-,1/2 

K  f»  fl*  X1J  Kx.x,  +  x-x.) 

/  /  e  rZ  dx.  dx, 

J-  J-uz-4iin 


r *  /•2*  i 

yo  yo 

-/*  YJ0<ST>  dT 


isY  cos (0-4) 


*  -j2.  J^(sa)  , 


Tfj  *  “wf  ^3  »  ^2*0. 


Following  the  same  procedure  as  In  Case  1,  expressions  are 
obtained  for  the  stress  components  of  the  forms 


-aP0  f*  y 2  cos  ♦ 


33  '  Jo 


[[(7  62  -  1)  hJ  cos2  ♦  -  rj]  wj  +  r3nj}  do  , 


J  ,  -aPo  ['  y2  s,n  ♦  ( 

12  —Jo 

. 


cosfc  ♦  w. 


1  r^3  ■  2y,y?)  7  7  .) 

2  - 2 — —  (sin*  $  -  cos*  ♦)  +  y3  W?  (  d#  , 

Y2  ) 


jf  .  ~aPo  /**  sin  2»  rt1  c'l 
23  ~TT  JQ  fl  YlY2tSl  “  S2-*  d*  * 


J31  ^T“  JQ  TT  ^ 


Y 2  COS*  +  Sx  -  (Y3  -  YjY2  sin*  ♦)  S2]  d*  ,  (66) 


where 


.  fm  -sy^x., 

St  •  /  J3(sa)  cos( sy  cos(e  -  ♦))  e  ds  , 

'  fm  -SYtX, 

Wt  «  /  JjC sa)  s1n(sY  cos(9  -  ♦))  e  z  *  ds  , 


t  »  1,2  . 

1  1 

The  Integrals  for  St  and  Wt  can  be  solved  analytically  so  that  the 
stress  expressions  have  the  form  of  single  integrals  which  can  be  eval¬ 
uated  numerically. 


Thermal  Stress;  Field 


Fi  and  72  are  given  by 


r  -  -W  [*  f 

1  J 


Ca2-x2]"2 


-[a2-,2]1'2 


cos(xjX^  +  XgXg)  dx2  dxj 


n*PnV  /*a  r 

**-* L  / 


r .2  2 ,1/2 

a  /**•  xl-*  i  Cx^x^+x^x,) 

/  e  1  L  £  £  dx,  dx, 

a  /-[a2  -  x2]1'2 


"MfPoV  aJi(sa) 
— E - s -  * 


2  -&r~ 


mt 


[a2-,2]1'2 

[a2-x2]2^2 


sinOT^Xj  +  x^)  dx2  dx^ 


J  sinfx^)  sin  (x2V?^xf  j  dx1 


Substitution  of  7i  Into  (61)  and  again  application  of  the  inverse 
transform  results  in 


°n  “2u*P.V  f •  /*- 

-  ■  -vrJo  Jo f2  cos(72X2) 


e  ■’(bjOj  +  b2D2) 


+  e  3(b,D,  +  b- D0)l  dx,  dx. 


The  expressions  of  the  stress  components  are  similar  to  those  (66) 

II  II  It 

in  Case  2  except  St  and  Wt  being  replaced  by  St  and  which  are  of  the 


s 


cos(sy  cos(Q  -  ♦))  e 


l  o 


ds 


* 


C  ' 


„  /*-  J?(sa)  -sYtx, 

Wt  *  4 y  -=-g —  sin(sY  cos(e  -  ♦))  e  J  ds  , 


t  ■  1,2  . 


It  II 

Again  St  and  W^.  can  be  solved  analytically  and  accordingly  numeri 
cal  results  of  the  stress  components  can  be  readily  obtained  by  auto¬ 
matic  integration  schemes  for  single  integrals. 


Thermal  Stress  Field 
It  is  found  that 


.[a2-x2l3l/2 


-UfP  V  f a  “X1J  /  y2\ 

1  «  l/2  i1  -  7)  +  ^  dX2  dXl 


’4wfPoV  Vsa) 
— R - ~T~ 


J2  *  0  . 

The  thermal  stress  field  can  be  obtained  by  replacing  f2  in  (67) 
with  f 3  which  has  the  following  form 

4J2(sa) 


As  an  Illustration  of  the  numerical  evaluation  of  the  Integrals, 
graphical  results  are  presented  in  dimensionless  form. 


The  coordinates  are  normalized  by  the  half  width  of  the  asperity  (a), 
such  that  (5,  n,  5)  =  (xi/a,  X2/a,  X3/a).  The  stress  components  are 
normalized  by  the  average  pressure  pQ,  0*^  *  o.j/pQ.  1116  temperature 
rise  from  the  cold  state  is  ♦  =  Tk/q0a,  where  q0  is  the  heat  flux 
due  to  the  average  frictional  load,  that  is,  q0  =  UfP0V.  Figures 
3-4  show  the  mechanical  principal  stresses  on  the  trailing  side  with 
different  depths,  ?,  for  Cases  1  and  2,  respectively.  It  is  noticed 
that  stresses  rise  sharply  In  the  neighborhood  of  S  »  1,  as  c 
approaches  zero.  Such  stress  singularities  at  5  3  ±1  are  expected 
owing  to  the  discontinuities  of  the  loading  conditions.  Temperature 
fields  for  different  cases  are  shown  in  Figures  5-7,  in  which  the  tem¬ 
perature  decays  rapidly  with  respect  to  the  depth.  Consequently,  the 
stress  component  on  is  compressive  only  in  a  very  thin  surface  layer 
due  to  the  thermal  effect.  Its  magnitude  of  compressive  stress  also 
decreases  rapidly  with  Increase  of  depth,  until  finally  it  becomes 
tensile.  The  consequential  effect  shown  in  Figure  8  that  the  maximum 
tensile  thermal  principal  stress  Is  at  c  »  10" 1  and  of  considerable 
magnitude.  It  should  be  noted  that  the  difference  between  the  magni¬ 
tude  of  the  thermal  stresses  (Fig.  8  here  and  that  In  [6])  is  due  to 
o^  In  [6]  being  the  dimensionless  normal  stress  Instead  of  the 
principal  stress.  The  combined  stress  field  is  represented  by  the 
maximum  principal  stress,  shown  in  Figure  9  for  several  depths  in  the 
near-surface  region.  Observe  the  the  Fig.  9  In  [6]  and  that  in  this 
paper  show  the  difference  only  of  order  2.  It  Is  attributed  to  the 
difference  of  2-0  and  3-0  analyses.  The  critical  asperity  pressure  is 
purposely  chosen  to  be  the  sane  (365  MPa).  To  Initiate  crack,  the 
asperity  travels  at  15  m/s  ("600  ips )  In  the  2-D  case,  while  the 
traverse  speed  has  only  to  be  10  m/s  ("400  ips)  in  the  3-0  case. 
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Trailing  Disiance  /  Radius  Of  Confacf  Area 


Figure  5.  Temperature  field  In  the  surface  layer  (Case  1). 


Figure  8.  Maximum  thermal  principal  stress  In  the  surface  layer  (Case  1). 


5.0  CONCLUSION 

The  paper  developes  the  general  mathematical  model  In  three* 
dimensional  formulation  to  simulate  the  thermomechanical  effect  near 
the  seal  surface  as  a  result  of  a  single  moving  asperity.  In  the  gen¬ 
eral  solution,  there  is  no  restriction  to  the  traverse  and  rubbing 
speeds  of  the  surface,  that  Is  the  traverse  speed  of  the  asperity  may 
not  be  tied  to  any  of  the  mating  surfaces  of  the  seal  assembly.  The 
contact  areas  by  postulation  of  the  present  analysis  have  the  shape  of 
rectangular  or  circular  regions.  However,  as  evidenced  from  the 
results  of  uniform  force  distributions  over  rectangular  and  circular 
areas,  the  shape  of  the  contact  area  is  not  that  critical.  With  refer¬ 
ence  to  Figure  5,  the  maximum  temperature  of  ♦  *  40  x  10” 3  corres¬ 
ponds  to  800*C  in  a  steel  based  seal  subjected  to  an  asperity  pressure 
of  365  MPa  (53,000  psi)  traversing  at  a  speed  of  10  m/s  (400  ips).  The 
same  temperature  may  of  course  be  attained  with  a  lower  asperity  pres¬ 
sure  of  240  MPa  (35,000  psi)  but  traversing  at  a  speed  of  15  m/s  (600 
psi).  The  high  temperature  is  maintained  only  near  the  surface  but 
drop  off  at  c  *  10“ For  the  data  used,  the  thermal  boundary  layer 
is  at  most  of  the  order  of  100  microns.  In  the  present  analysis,  it  is 
postulated  that  the  mechanical  properties  of  the  seal  material  remain 
unchanged  for  all  temperature  levels  during  the  loading.  For  those 
seal  materials  that  m^y  change  their  mechanical  properties  at  the  high 
temperatures  near  the  boundary,  the  boundary  layer  analysis  needs  to  be 
considered.  It  is  not  unexpected  that  the  maximum  tensile  stress 
should  occur  at  the  thermal  layer  Interface.  The  subsurface  initiation 
of  fracture  is  also  found  In  the  two-dimensional  model  at  a  depth  of 
one-twentieth  of  the  asperity  width  [7].  The  average  friction  required 
for  fracture  initiation  is  still  high.  It  corresponds  to  a  pressure 
around  240  MPa  (35,000  psi)  in  Figure  9,  and  a  high  Coulomb  coefficient 
of  friction  Uf  *  0.5. 
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